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Abstract 

We construct a Schrodinger-like equation for the longitudinal wave function of a meson in the 
valence qq sector, based on the 't Hooft model for large-iV two-dimensional QCD, and combine 
this with the usual transverse equation from light-front holographic QCD, to obtain a model for 
mesons with massive quarks. The computed wave functions are compared with the wave function 
ansatz of Brodsky and De Teramond and used to compute decay constants and parton distribution 
functions. The basis functions used to solve the longitudinal equation may be useful for more 
general calculations of meson states in QCD. 



PACS numbers: 11.25.Tq,12.38.Lg,12.39.Ki,14.40.-n 



I. INTRODUCTION 



Light-front holographic QCD [l| provides a Schrodinger-like equation for the transverse 
wave functions of hadrons, based on a hypothesized duality with AdS 5 . The duality argu- 
ments lead to the construction of models for a confining transverse potential. The soft-wall 
model [2| in particular admits analytic solutions for the spectrum and the transverse wave 
functions, and, what is more, yields a spectrum 0-0] consistent with linear Regge trajecto- 
ries. However, because the duality relies on the approximate conformal limit of zero-mass 
quarks, the transverse models do not include dependence on quark masses. Also, the lon- 
gitudinal wave function is only specified for zero-mass quarks, through a matching of dual 
form factor calculations 

One escape from these restrictions is an ansatz by Brodsky and De Teramond 0] that 
extends the transverse momentum dependence to include the full invariant mass. For meson 
states, this extends an exponentiation of kj_/x(l — x), where k± is the transverse momentum 
and x is the longitudinal momentum fraction, to k\/x(l — x) + m\/x + m\j{\ — x). To 
be consistent with the zero- mass limit, we take the quark masses mi and 1712 to be current 
quark masses. 

An alternative is to assume separation of variables for the meson wave function and then 
write down a light-front equation for the longitudinal part that includes quark masses and 
a model potential. The kinetic term will be just m\jx + m\j{\ — x), though here, to be 
consistent with nonrelativistic quark models, the quark masses are constituent masses. The 
potential term should be confining and should yield longitudinal wave functions consistent 
with the dual form-factor analysis, in the zero-current-mass limit. 

A potential model that achieves this, and is directly related to QCD, is the 't Hooft 
model obtained in the large- A limit of two-dimensional QCD . There exist nearly exact 
analytic solutions for the ground state, which can be improved easily with numerical calcu- 
lations (ll - lp. In fact, a natural basis set for the numerical calculation, which uses Jacobi 



polynomials [13(, could be useful for general calculations of meson valence states. The cou- 
pling constant in the model is constrained by the requirement that the longitudinal wave 
function for the pion be consistent with the dual form-factor analysis when the constituent 
quark masses have appropriate values. 

Here we explore this potential-model approach. The details of the model and its solution 
are discussed in Sec. [Hi Sample applications are illustrated in Sec. IHH A summary and some 
additional remarks are included in Sec. [TV] 



II. THE MODEL 

In light-front holographic QCD, the meson valence wave function is written in factorized 
form 

ij nLS (x, C, <p) = X(x)cf> nLS (0-^=e iL v, (2.1) 



where n is the transverse radial quantum number, L is the orbital angular momentum, S is 
the spin, and £ = x(l — x)b± is the transverse variable, with b± the transverse position. 
The transverse wave function (f> n LS satisfies a Schrodinger-like equation, which in the soft- 
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wall model, takes the form [l|, 0] 
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with Mi? the longitudinal-transverse separation constant. The analytic solutions are two- 
dimensional oscillator wave functions 
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where the L% are associated Laguerre polynomials. The mass spectrum is 
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In our approach, the longitudinal wave function X(x) satisfies its own Schrodinger-like 
equation 
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with T 7 indicating the principal value and C a constant. This is the equation obtained in 
the 't Hooft model for large- iV two-dimensional QCD lb] ]. Because the transverse equa- 
tion already introduces enough quantum numbers, we consider only the ground state of the 
longitudinal equation; any additional quantum number associated with longitudinal exci- 
tations would represent double counting. Also, since the light-meson spectrum is already 
represented by the transverse equation, the constant C is used to set My to zero. The net 
effect is that the additional longitudinal equation is only for determination of the longitu- 
dinal wave function and has nothing to say about the spectrum. To recover nonrelativistic 
quark models in the appropriate limit, we use constituent masses for m\ and m 2 . 



As discussed by 't Hooft [10 



and Bergknoff lH, the ground-state wave function X(x) is 
well approximated by the form x^(l — %Y 2 . Their analyses of the endpoint behavior show 
that Pi should satisfy the transcendental equation 
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From the dual form-factor analysis [8] , the light-meson wave function should have this form 
with 0i = 02 = 1/2. If we take the up and down quark masses, m u and mj, to be equal, the 
square-root behavior is obtained if g 2 /7t = m 2 , consistent with cot7r/2 = 0. This fixes the 
value of the coupling constant. Although the model could be more flexible if g were flavor 
dependent, we do not consider this. 

The exact solution cannot be obtained analytically. Instead, we solve (12.51) numerically 
by expanding the wave function X in terms of orthonormal basis functions /„, chosen to 
include the analytic approximation explicitly. As discussed by Mo and Perry [l3||, these 
basis functions are 

f n (x) = N n x^(l - x)^P^ 2h \2x - 1), (2.7) 
with Pn 2 ^ 2 ' 213 ^ the Jacobi polynomial of order n. The normalization factor N n is given by [1 
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The solution is then represented as 



X(x) = Y,Cnfn(x). (2.9) 

n 

For equal- mass cases, the longitudinal equation obeys an x -H- (1 — x) symmetry, and only 
the even-n terms will contribute. In general, we find that only a few terms are needed; the 
n = term, for which the Jacobi polynomial is constant and fo oc x^(l — x)^ 2 , represents 
90% or more of the probability. 

The expansion coefficients c n are obtained by diagonalizing the longitudinal equation in 
the /„ basis. The matrix representation is 

(TH/^ 777. ^ \ 

-^A X + -%A 2 + B c = £c, (2.10) 

with £ = C/mi and 

f 1 dx f 1 dx 

(Al) nm = / f n (x)f m (x), {A 2 ) n m = / Z fn{x)fm{x), (2.11) 

Jo x Jo 
B nm = f dxvf ' dyf n {x) Sm{ f~ S ^ V \ (2.12) 

Jo Jo \ x ~ y) 



Following 't Hooft [lOj, the matrix representation is made explicitly symmetric by rewriting 
the potential term as 

2 Jo Jo x-y x-y 

In addition to the explicit symmetry, which simplifies the matrix diagonalization, this rear- 
rangement also resolves the principal value prescription. The matrices are small because the 
number of terms needed in the expansion are few; the diagonalization is then straightforward. 



III. SAMPLE CALCULATIONS 

To see the implications of our model, we compare the form of the longitudinal wave 
function with the ansatz by Brodsky and De Teramond Q , both directly and through the 
computation of parton distribution functions and decay constants, for the pion, kaon, and 
J/ty. Where parameter values are needed, we use the current-quark parameterization of 
Vega et al. jl| with no additional fits or adjustments. The parameter values are listed in 
Table E 

For the ground state, the complete wave function is 

ij(x,() = NX(x)e- R2<:2/2 (3.1) 
or, in terms of the transverse coordinate bj_, 

^(ar, & l) = NX(x)e- K2x{1 - x)b ^ /2 . (3.2) 
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TABLE I. Meson parameters and decay constants. All dimensionful parameters are in units of 
MeV. Results are compared between our model, which uses constituent-quark masses, and the 
ansatz of Brodsky and De Teramond 

experimental values are from Vega et al. 4j and the Particle Data Group (l5j. 



, where we use current-quark masses. Parameter and 



model ansatz decay constant 

meson mi ui2 mi 7712 P qq K model ansatz exper. 
pion 330 330 4 4 0.204 951 131 132 130 
kaon 330 500 4 101 1 524 160 162 156 
J/tf 1500 1500 1270 1270 1 894 267 238 278 



The factor X is fixed by the normalization 
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where P qq is the probability of the quark-antiquark valence state. If X(x) is separately 
normalized such that 

^« = l, (3.4) 



(3.5) 



x(l — X) 

then X = - \/P qq - A Fourier transform to transverse momentum k±_ yields 



ip(x,k±) 
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For zero-mass quarks, X(x) = \fxji — x), in order to satisfy the form factor duality [8(. 
The ansatz for X(x) by Brodsky and De Teramond [9[ extends this to massive quarks by 
allowing k\/x(l — x) to become the invariant mass {m\ + k*j_)/x + (m| + k\)/(l — x) with 
the addition of the two mass terms. This gives 



with 
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Figure [T] compares this ansatz with the X(x) computed in our model. For the pion, the 
two wave functions are essentially the same, since both involve only tiny variations from 
y/x{l-x). 

The wave functions can be used to compute decay constants, as given by [1( 



f M = 2^6 [ dx f 
Jo Jo 



dk\ 



(3.8) 



The results are included in Table U The values for the pion and kaon are consistent and in 
agreement with experiment. The value for the J/\l/ is significantly closer to experiment than 
the value obtained from the longitudinal ansatz. 
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FIG. 1. Longitudinal wave functions X{x) for the (a) pion, (b) kaon, and (c) J/^. The solid lines 
are wave functions from our model; the dashed lines show the ansatz by Brodsky and De Tera- 



We can also compute parton distributions f(x). As discussed in Q and shown in 17 . if 
the wave function takes the form 

1>(x, k ± ) = r/(x)e- fe i /(2K2;c(1 - x)) , (3.9) 

the parton distribution is given by 

f(x) = -^x(l-x)r ] 2 (x). (3.10) 



Applying this to our model, we obtain 

•\ = P- 

x(l — x) 



m = p m^K- (3-11) 
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For the ansatz, we have 

W*) = Nl^P qq e-^* +m y^l«* '. (3.12) 

These are plotted for comparison in Figure [21 The rough similarity of the wave functions 
translates into similar parton distributions. 

0.4 n 1 



0.3 - 




(b) (c) 
FIG. 2. Same as Fig. CD but for parton distributions f(x) multiplied by x. 

IV. CONCLUDING REMARKS 

We have constructed and solved a relativistic light-front equation for the longitudinal 
wave functions of mesons with massive quarks, to be used in tandem with the transverse 
equation of light-front holographic QCD. Comparisons with the ansatz [9( (13. 6p show that 
for lighter mesons, the longitudinal wave functions are quite similar. However, for the J/ l I / , 
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there is a notable difference, which translates into a better estimate of the decay constant, 
as listed in Table HI Wave functions for the J/\l/, and also the pion and kaon, are shown in 
Fig. [TJ and parton distributions are shown in Fig. [2j The similarity of the longitudinal wave 
functions provides the ansatz with a connection to the fundamental interactions of QCD. 

Perhaps the most broadly useful outcome of this exercise is to illustrate that there is a 
convenient set of basis functions for longitudinal wave functions of meson valence states. 
These are the f n defined in (12. 7p . with the parameters /?, to be optimized as needed for a 
given application. One could, of course, use the eigenfunctions of the 't Hooft model, but 
these do not have an analytic form and would add an extra layer of complication to any 
calculation. 

The importance of the choice of basis functions is in the rate of convergence as the 
basis is expanded. For the alternative of discrete light-cone quantization, it is known that 
convergence can be much slower [l8|. Thus, these basis functions may prove useful for 



calculations, such as those described in [19|, that use the transverse light-front holographic 
eigenfunctions. 
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